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1 Introduction 

Some important models of quantum field theory and condensed matter physics exhibit 
remarkable properties granting them the gift of integrability. Among prominent examples 
we have the quantum non-linear Schrodinger model [I] , the Heisenberg spin chain [2 j and 
the Hubbard model j3). The realm of those models are of quantum nature and, although 
a rigorous and unambiguous definition of integrability at quantum level is still lacking 
(4^[T] , these models exhibit a set of enhanced symmetries which can be explored in order 
to compute physical quantities exactly Among those quantities we mention the 

energy spectrum and some correlation functions, and in this way integrability offers non- 
perturbative access to the behaviour of strongly interacting systems. These achievements 



are in large part due to the advent of the Bethe ansatz 10 and its algebraic formulation 



1112 characterising the model wave function in terms of its momentum. This algebraic 



formulation have been successfully applied to the majority of known integrable systems 



13 -16 , yielding the model eigenvectors in terms of creation operators acting on a suitable 
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reference state. Moreover, as remarked in [17] the study of Bethe's wave function also 



provided important insights leading to the concept of commuting transfer matrices 18 



On the other hand, the exact solution of a model whose transfer matrix belongs to a 



commutative family is also intimately connected with functional equations methods 19 



As far as spectral properties are concerned, we have available a variety of functional 
methods |19)-22 yielding transfer matrices eigenvalues which do not address the problem 
of explicitly constructing the respective eigenvectors. In fact, the evaluation of most 
physical quantities do not require the eigenvectors themselves but quantities which can 
be derived from them. For instance, the calculation of correlation functions would then 
require the evaluation of eigenvectors scalar products in addition to the expected value 
of operators [8]. 

Within the framework of the algebraic Bethe ansatz, the calculation of scalar prod- 
uct of Bethe vectors was inaugurated by Korepin 23 largely influenced by Gaudin's 



hypothesis on the norms of the non-linear Schrodinger model wave function 24 . The 



method of 23 has also been described in [8] where the scalar product of Bethe vectors for 
the six-vertex model is given as the vacuum mean value of a determinant whose entries 
are expressed in terms of quantum fields. In [8] that scalar product is also given as a 
summation over the product of determinants with scalar entries. 

The purpose of this paper is to demonstrate that scalar product of Bethe vectors 
can also be computed by means of functional equations. The method we shall employ 
here resembles the one described in the series of works 25 -28 where the Yang-Baxter 
algebra has played the major role in deriving functional equations describing the partition 
functions of vertex and SOS models with domain wall boundaries. The origin of the 
functional equations describing scalar products is also the Yang-Baxter algebra and its 
solution turns out to be given by a multiple contour integral. 

This paper is organised as follows. In Section [2] we give a brief description of Bethe 
vectors for the six- vertex model. In Section [3] we illustrate how the Yang-Baxter algebra 
can be explored in order to derive functional equations describing scalar products. The 
solution of our functional equations is obtained in Section [4] and concluding remarks are 
discussed in Section [5| Technical details and proofs are presented through the Appendices 
A to E. 



2 Bethe vectors for the six-vertex model 



Vertex models in two dimensions constitute one of the corner stones of the theory of 
exactly solvable models of statistical mechanics 2! 
eigenvectors of the six-vertex model |30] by means of the Bethe ansatz 



In particular, the study of the 

provided 

insightful information which paved the way to establish the connection between two- 
dimensional vertex models and one- dimensional spin chains 31 , 32 , 18 . This study also 



received a large impulse with the advent of the Quantum Inverse Scattering Method 
(QISM) which unveiled the algebraic foundation supporting the construction of Bethe 
vectors [TT , 12 . 
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Bethe vectors. In the framework of the QISM, the eigenvectors of the six-vertex 
model are built up from the action of creation operators on a pseudo-vacuum state. 
More precisely, an eigenvector is of the form 

\if>) = B{\»)...B(\»)\0) (2.1) 

where B(X), together with three more generators A(X), C(X) and D(X), satisfy a certain 
set of algebraic relations. In its turn, the vector |0) is the sl 2 highest weight vector while 
the parameters A^ G C are required to satisfy certain constraints. 



Yang-Baxter algebra. The algebra satisfied by the generators A(X), B(X), C(X) and 
D(X) is commonly referred to as Yang-Baxter algebra and it reads 



fti 2 (Ai - A 2 )Ti(A 1 )r 2 (A 2 ) = r 2 (A 2 )Ti(Ai)^i 2 (A 1 - A 2 ) 



(2.2) 



The relation g2) is defined in End(Vi <g> V 2 ) with Yj = C 2 . In their turn, Ti = T ® id 
and % = id £g> T with 



(A(X) B(X)\ 
' [A) \C{X) D(X) J 



(2.3) 



The matrix TZ a b G End(V a <8> V&) encodes the algebra structure constants associated with 
the six-vertex model and it is explicitly given by 



n 



fa 0\ 

b c 

c b 

\0 a) 



(2.4) 



where a(A) = sinh(A + 7), 6(A) = sinh (A) and c(A) = sinh(7). The 7£-matrix (2.4) 
satisfies the Yang-Baxter equation 

ftia(Ai - A 2 )^i 3 (Ai - A 3 )^ 23 (A 2 - A 3 ) = ft 23 (A 2 - A 3 )^ 13 (A 1 - A 3 )^i 2 (Ai - A 2 ) (2.5) 



in End(Vi <8> V 2 <8> V 3 ) ensuring the associativity of (2.2). 



Representations. The representations of A(X), B(X), C(X) and -D(A) in the tensor 
product space Vi ® • • • ® Vl can be built as follows. We consider an ordered product of 
^-matrices in the tensor product space V a ® Vi <E> • • • <8> Vl, namely 



f a {x)= J] n aj {x- N ) 

1<0<L 



(2.6) 



with parameters X,fij G C. Due to (2.5) one can show thatT a satisfy (2.2) and thus 
T =T a yields a representation of A(X), B(X), C(X) and D(X). 



4 



Highest weight vector. The vector |0) in (2.1) is the 5^ highest weight vector and 
it explicitly reads 



10) 



(2.7) 



in the space Vi 



Yl- Taking into consideration the representation (2.6) with 



7?.-matrix given by (2.4), we readily obtain the properties 

L 



A(X)\0) = l[a(X-^)\0) 



A(\) = l[a(\- N )(0\ 



3=1 
L 



3=1 
L 



D(X)\0) = l[b(X-f, J ) |0> 

i=i 

C(A) |0) = 

where (0| corresponds to the transposition of |0). 



(0\D(X) = l[b(X- f i j )(0\ 

3=1 

(0|B(A) = 0, 



(2- 



Twisted transfer matrix. The Yang-Baxter algebra (2.2 ) enables us to show that the 



matrix T(A) = (j)iA(X) + 02-D(A) forms an one-parameter family of mutually commuting 
matrices, i.e. [T(Ai), T(X2)} = 0. The matrix T(X) is called twisted transfer matrix and 
for vertex models it represents the configuration of rows spanning a two-dimensional 
lattice with periodic boundary conditions. In their turn, the parameters 0i,02 £ C 
govern the deviations from strict toroidal boundary conditions and, in this sense, they 
introduce the notion of twisted boundary conditions |33|. The vector \ip) as defined by 



(2.1) is an eigenvector of T(A) with eigenvalue 



Aw^^naix-^n 

i=l ^ 1 



o(Af - A) 



i=i 



A) 



+ 



q(A - Af ) 
fi & (A - Af ) 



(2.9) 



for particular choices of the parameters Af 



Dual Bethe vector. The dual vector 

M = <0|C(Af)...<7(A?) (2.10) 

is an eigenvector of T, i.e. T(A) = A(X) with the same eigenvalue A as given by 
(2.9). Nevertheless, here we shall consider the dual vector 

£ = <0| C{\? ) . . . C{X c n ) (2.11) 



in order to keep our results as general as possible. 
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Off-shell scalar product. The scalar product of Bethe vectors S n is then defined as 

n n 

S n (\?,...,\Z\\?,...,\*) = (0\l[C(\f) n^(\ B )|0) • (2-12) 



i=i 



We shall refer to (2.12) as off-shell scalar product when the variables Af and Af are free 



to assume any value on the complex plane. 



On-shell scalar product. The vector defined in (2.1) is an eigenvector of the 
twisted transfer matrix T(A) only when the variables Af satisfy certain constraints. 



Thus the parameters Af and Af need to be fine tuned in order to having (2.12) de- 
scribing the norm of a Bethe vector. The constraints on the variables Af ensuring that 



(2.1) is a transfer matrix eigenvector are the so called Bethe ansatz equations. Similar 



constraints would be required for the dual eigenvector (2.11) but here we shall consider 
a less restrictive condition by keeping the variables Xf arbitrary. In this way the scalar 



product (2.12) under the condition 



n 



-1<P2 



«(Af - Af ) 



(2.13) 



k 



will be referred to as on-shell scalar product. 



3 Functional equations 

This section is devoted to the derivation of functional equations describing the scalar 



product S n defined in (2.12). The main ingredient of our derivation is the commutation 
rules encoded in the relation (2.2) commonly referred to as Yang-Baxter algebra. We 



shall obtain two different functional equations and the determination of S n will rely on 
the resolution of this system of equations. This method consists of an extension of the one 



originally proposed in 25 for the partition function of the six-vertex model with domain 
wall boundaries and subsequently generalised in 27 28 for SOS models. Similarly to 



those cases, the derivation of functional equations for scalar products will explore a 
consistency relation between the si2 algebra highest weight representation theory and 
the Yang-Baxter algebra. 

3.1 Equation type A 

We consider the quantity 



(0\i[C(Xf)A(X )l[B(Xf)\0) 



(3.1) 



i=l 
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computed in two different ways. In the first way we consider the commutation rules 
B(X 1 )B(X 2 ) = B{X 2 )B{\ 1 ) 



A(X 1 )B(X 2 



a(X 2 - Ai) , . , c(A 2 - Ai) 



b(X 2 - X 1 ) 



B(X 2 )A(X 1 ) - 



b(X 2 -X 1 ] 



B(X 1 )A(X 2 ) 



(3.2) 



contained in the relation (2.2) to move the operator A(Xo) in (3.1) to the right through 
the string of operators B(X~f). Due to ( |3.2 ), at the last step we shall need to compute 



the action of a given operator A(X) on the state |0). For that we use the relation (2.8) 
arising from the st 2 highest weight representation theory. By doing so we find that the 



quantity (3.1) consists of a linear combination of terms S n with 5-type arguments in the 
set {A , Af , . . . , Af } with only n elements being taken at a time. 

The second way of evaluating (3.1 ) is by moving the operator A(Xq) to the left through 
all the operators C(A^ 



C(X 1 )C{X 2 ) 
C(X 1 )A{X 2 ) 



This can be implemented with the help of the relations 
= C(X 2 )C(Xx) 



a(Ai - A 2 ) 
6(Ai - A 2 ) 



A(X 2 )C(X 1 ) - 



c{Xi - Aj 
6(Ai - A 2 ) 



A(Ai)C(A 2 



(3-3) 



which are also among the ones encoded in (2.2). Then at the last stage we will need the 



quantity (0| A(X) which is given in (2.8). Thus this alternative route of computing (3.1) 



also yields a linear combination of terms S n but this time with C-type arguments in the 
set {A , Af , . . . , Af} where only n variables are taken at a time. 



In this way the consistency between these two routes of computing (3.1) implies the 
functional equation 



M S n (Xi , . . . , Af |Af , . . . , Af ) + N- S n (Xi, . . . , A^|Aq, Af , 



A i+li 



An) 



+ ^2 S n (Xo, Af , . . . , Xf_ x , Xf +1 , . . . , Af|Af , . . . , Af) — 



t=i 



with coefficients 



N 



(B,C) 



c(Af' c 



n 



a(Af - A ) 



\ HAf - A ) t\ b(Xf - A ) 



n 



a(Af - A ) 



. x L n (\B,C \B,C\ 

AoJ TT (X B,C \ TT a \*j ~ A i ) 



KK - Ac 



3=1 



6(A, - X { ) 



where as = 1 and «c = — 1. 



(3.4) 



(3.5) 



3.2 Equation type D 



The same mechanism employed in Section 3.1 can also be considered starting with the 
quantity 



(0\HC(Xf)D(X o )l[B(Xf) |0> 



(3.6) 



i=i 
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instead of (3.1 ). In that case, moving the operator D(\ ) to the right through the string 



of operators B(Xf) will require the use of the Yang-Baxter algebra relation 
D(X 1 )B(X 2 ) = ?^ " ^ BMDfr) - ~ X "\ b{X 1 )D{X 2 ) 



b(Xi - A 2 ) 

On the other hand, the commutation rule 

a(A 2 — Ai 



6(Ax-A 2 ) 



(3.7) 



C(X l )D(X 2 ) 



b(X 2 - Ax) 



^(A 2 )^(A 1 ) - Xl) D(X 1 )C(X 2 ) 
b{X 2 - Ai) 



(3i 



will be required in order to move the operator D(Xq) to the left through all the operators 
C(Xf). Besides that we shall also consider the commutation rules [B(Xi), B(X 2 )] = 



[C(X 1 ),C(X 2 )] = and the properties (2.8). 



Thus the consistency condition between these two different ways of computing (3.6) 
leave us with the following functional equation, 



M S n (Xi , . . . , A^|Af , . . . , Af ) + N- S n (Xi, . . . , A^|Aq, Af , 



i=l 



+ ^ Ni° S n (Xo, Af , . . . , Xf_ x , Xf +1 , . . . , A^|Af , . . . , Af ) — 



i=l 



whose coefficients are explicitly given by 

L 



Mi 



o 



(B,C) 



Y[b(X -fij) 

3=1 



A i-li A i+li 



tt a(Ao - Af) A a(A - Af ) 
ll6(Ao-Af) f = {6(Ao-Af) 

B,C \B,C\ 



OCB,C 



/\ \B,C\ L n x 



Af ) 



6(A - Af - c ) f* 



, < &(Af '° - Af ' C ) 



(3-9) 



(3.10) 



In summary, we have demonstrated in this section how the Yang-Baxter algebra 
can be explored in order to derive functional equations for the scalar product of Bethe 
vectors. We have obtained two distinct equations which we shall refer to as equation 



of type A (3.4) and equation of type D (3.9). The solution of these equations will be 



discussed in the next section. 



4 Scalar product 



Solving the system of functional equations formed by (3.4) and (3.9) is the main goal 
of this section and some remarks are required in order to proceed. For instance, the 



method employed here for the derivation of (3.4) and (3.9) can be seen as an extension 
of the method considered in 28 , and the resulting functional relations indeed share 



some similarities with the one obtained for the partition function of the elliptic SOS 
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model with domain wall boundaries. However, there are still some important structural 



differences that introduce some extra steps in solving (3.4) and (3.9) 



Firstly, let us consider the similarities. The functional equations obtained here and 
the one of 



28 are all relations for a multivariate function, i.e. F{ 



,z n ), composed 

of a linear combination of terms involving the function F with a given variable z% in its 



argument being replaced by a variable z . Moreover, the Eqs. (3.4) and (3.9) are also 



homogeneous in the sense that aS n solves our system of equations if S n is a solution and 
a is a constant. In fact, a only needs to be independent of the variables Af and Af . This 



property tells us in advance that Eqs. (3.4) and (3.9) will be able to determine S n only 
up to an overall multiplicative factor independent of Af and Af at most. Thus it will 
be necessary to evaluate the function S n for a particular value of its variables in order 
to have our scalar product completely fixed. Still considering the similarities, both Eqs. 
(3.4) and (3.9) are linear which raise the issue of uniqueness of the solution. Here we will 



be interested in a multivariate polynomial solution and this property ensures uniqueness 
as demonstrated in |27] under very general conditions. 



Now let us consider the differences between (|3.4| 3.9|) and the functional equation 
obtained in 



2S 



The most obvious difference is that here we have obtained two equations 
which might suggest that one of them is redundant. However, the direct inspection of 
our equations for small values of n and L reveals that the polynomial solution which 
we shall be interested can not be completely fixed by only one of the equations. The 
situation is different when we consider both equations simultaneously, and their direct 



inspection shows that the system (3.4, 3.9) indeed determines the polynomial solution 

in the case considered in 28 we 
Ai, . . . , A n }. Here both 



up to an overall multiplicative factor. Furthermore 
have an equation running only over one set of variables, i.e. {Ao 
of our equations run over the two sets of variables {Ao, Af , . . . , Af } and {Ao, Af , . . . , Af }. 
Taking into account the above discussion, the following Lemmas will pave the way for 



solving (3.4, 3.9). 



Lemma 1 (Polynomial structure). In terms of variables x 



B,C 



[x 



is of the form S n = Yl^ =1 (xfxf) 2 S, 
of order L — 1 in each one of its variables separately. 



o 
1 • 



1 ) 



X 



B\ 



Proof. See Appendix [A] 

Lemma 2 (Special zeroes). The function S^Af , 
specialisation of variables Af = jii and Af = /ii — 
the specialisation Af 



= e 2Xi ' , the function S n 
where S n is a polynomial 

□ 



fii and Af 



. . . , Af |Af , . . . , Af ) vanishes for the 
7. The same property also holds for 



Mi -7- 



Proof. See Appendix I 



□ 



Lemma 3 (Doubly symmetric function). The scalar product 5* n (Af , . . . , Af |Af , . . . , Af ) 

is a symmetric function in each one of the set of variables {Af } and {Af } independently. 
More precisely, 



S n ( A i 



c 



A 



c 



A n \ A l J 



— S n (X 1 , 



Af |Af , . . . , Af ) 
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and 



S n (X l , . . . , X n \X 1 , . . . , Aj , . . . , Xj , . . . , A n ) 



C (\C \V\ ^B x b )Hn 

^n^l , • • • , * n \ A 1 , . . . , Aj , . . . , Aj , . . . , A n ) . 



,C| \B 



Proof. See Appendix O 



□ 



Lemma 4 (Asymptotic behaviour). In the limit xf ,c — > oo, the function S n behaves as 



l<ai<---<a n <L 

where [n\] q 2 denotes the q-factorial function defined as 

[n\} q2 = 1(1 + g 2 )(l + q 2 + q 4 ) . . . (1 + q 2 + ■ ■ ■ + g 2 ^ 1 )) . 



(4.1) 



Proof. See Appendix D 



□ 



Remark 1. Due to the Lemma\3\we can safely employ the notation 

S n (Xx, . . . , A^|Af , . . . , Af ) = S , n (X 1,?1 |F 1 ' ?1 ) 
where X iJ = {X% | % < k < j} and Y l ' j = {Af \i <k < j}. 



Except for the Lemmas [T] and 



functional relations (3.4) and (3.9). 



4j the remaining ones are a direct consequence of the 



4.1 Off-shell formula 



Here we shall consider the resolution of the system of equations (3.4, 3.9) for general 



values of variables Af and Af. This case is refereed to as off-shell scalar product as 



remarked in Section 2 The methodology we shall employ for solving (3.4, 3.9) is similar 



to the one developed in 28 and in what follows we describe a sequence of steps leading 
to the solution. 



Step 1. We firstly consider Eq. (3.4) under the specialisation of variables A = /ii — 7 



and Af = fi± such that the coefficient M vanishes. We shall also consider the property 
S n (X 1,n \ni — 7, . . . = obtained from Lemmas and M implying the relation 



S n {X 2 ' n \Y 2 ^ 



n&(A 

J=2 



c 



Vi)a(Xf -/iO V(X 2 ' n \Y 2 ' n ) 



(4.2) 



due to Lemmas 1 and 3 In (4.2) we have also introduced the notation Z M = Z 1 ' 3 U{/Xi 



7} and Z 1 ' 3 = Z 1 ' 3 U {^1} for Z 1 ' 3 G {X l,3 \ Y 1 ' 3 }. In its turn the function V appearing 
in (4.2) is also of the form described in Lemma [l] under the identifications S n 1— > V, 
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ti i — y Ti — 1 and L i— > L — 1. Thus under this specialisation of variables, Eq. (3.4) yields 
the relation 



(4.3) 



i=i 



where X^ = X*>>\{\<j} and 



n-l 



m 



(A) 



N, 



(C) 



J=l 3=1 iV ™ 



n 



«(Af - /^) 



(Af - 0i) jLi a 0i - Mi) 7 = i 



3=1 



g(Af - Af ) 

KAf - Af ) 



(4.4) 



Step 2. Analogously to (4.2), due to Lemmas [1] |2| and |3j we can write 

n 

S n (X 2 ' n \Y 2 ' n ) = 11 a(Xf - /ii)6(Af - m) W(X 2 > n \Y 2 > n ) 

3=2 



(4.5) 



where the function W is also of the form described in Lemma [T] under the mappings 
ti i — y ti — 1 and L \-¥ L — 1. Then by setting A = H\ — 7 and A^ 7 = \i\ in Eq. (3.4), and 



considering the relation (4.5) in addition to the property S n (fii — 7, ... , fii\Y 1,n ) = 0, we 
obtain the formula 



S n {X 1 > n - 1 \Y 1 > n ) = ^m 4 (A) W(X l > n - l \Yl> n ) 
where Y^ = F ij \{Af } and 



(4.6) 



i=i 



n-l 



m 



(A) 



N. 



(B) 



n 



3=1 



A 0=M1~7 



a(Af - ^ 



n-l 



a (Af - Mi) y a(/xi - Hj) fJi 



n n(\ B — \ B \ 



HA? - Af ) 



(4.7) 



Step 3. The coefficient vanishes under the specialisation Af = /ii — 7. Thus, 

under this particular specialisation, the Eq. ( 3.4[ ) will contain only terms of the form 
S n (X 1,n \Y 1,n ~ 1 ) allowing us to use the relation (4.3 ) to obtain an equation involving solely 
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the function V. The resulting equation can be further simplified by setting A^ = fii, 
and in this way we are left with the relation 



n-l 



n-l 



Jo V{X l ' n - 1 \Y 1 ' n ~ l ) + J2 K i B)v ( xl ' n ~ 1 \ Y z°' n ~ 1 ) + ^Kf ) V(X° i ' n - l \Y 1 ' n - 1 ) = . 



i=l 



i=l 



(4i 



In their turn the coefficients Jo and K\ a ppearing in (4.8) correspond respectively to 

(B C) 

the coefficients M and iV> ' given in (3.5) under the mappings L \— > L — 1, n \— y n — 1 



and fii i — y /ij+i- Thus the function V obeys essentially the same equation as S n _i for 
Bethe vectors living in Vi <g) • • • <g) Vl_i. 



Step 4. Next we consider Eq. (3.4) with X% = Hi — 7. In that case the coefficient 
N n C ^ vanishes and we are left only with terms of the form S , n (X 1,ri ~ 1 |y 1,?1 ). Then we use 



the relation (4.6) to obtain an equation only in terms of the function W. After setting 



X B = Hi and eliminating an overall factor, the equation obtained in this way reads 



n— 1 n— 1 

*( B hxr( vi>n-liv°. n -l\ 1 jA c )\ 



Jo W(X 1 ' nr - 1 \Y 1 ' n - 1 ) + ^K\ B) W{X 1 > n - l \Y?> n - x ) + J2 K i W ( X i' n ~ 1 \ Yl,n ~ l ) = • 



i=l 



(4.9) 



Thus the functions V and W obey the same equation and the uniqueness of the solution 



of (3.4) would imply that V and W can differ only by an overall constant factor. 



Step 5. We consider the Eq. (3.9) with Ap = Hi an d X B = Hi — 7. In that case the 



coefficient M = 0, and we can use formula (|4.5|) in addition to the Lemma |2| to obtain 
the relation 



(4.10) 



i=l 



where 



ft-i 



- (D) 



N, 



(C) 



i=i 



A =Ml 
A n=A»l-7 



n 



Hp, - Af ) 



n-l 



]Ja(\f - ni)Y[a(\f - hi 



q(Af - xf) 



(4.11) 
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Step 6. Next we set A = \i\ and A^ = ^1—7 in the Eq. (|3.9[). This allows us to employ 



the relation (4.2) and the property S n (fii, . . . , /^i — j\Y ,n ) = described in Lemma \2 



In this way we obtain the expression 

n 

S n {X 1 ' n - x \Y 1 ' n ) = Y,mf ] V{X x ' n - l \Y X; 



(4.12) 



i=l 



with 



n-l 



m 



(£>) 



-n 6 (Af-^)n°( A f-^ 



iV, 



(B) 



3=1 



A', 



(C) 



n 



n-l 



6(^1 - A f ) f = i a (^' - Mi) } = i 



IJa(Af-Mi)II a ( A ?-^) 



«(Af - Af ) 



3 = 1 



6(Af - Af ) 



(4.13) 



Step 7. The Eq. (3.9) will contain only terms of the form S n (X ' n \Y ' n ) for the 



specialisation X? = \L\ since the coefficient vanishes. In that case we can consider 
the formula (4.10) and also set A^ 7 = — 7. By doing so we are left with the relation, 

n— 1 n— 1 

Jo W(X 1 ' nr - 1 \Y 1 ' n - 1 ) + Y,K i i B) W(X 1 ' n - 1 \Y^ n ' 1 ) + Y,W ) W{X°> n ~ 1 \Y 1 > n - 1 ) = , 

(4.14) 



i=i 



i=l 



whose coefficients Jq and iC ' correspond respectively to the coefficients M and N 



(B,C) 



given in (3.10) under the mappings L \- > L — 1, n h > n — 1 and /ij+i. 



Step 8. Set A^ = \i\ in Eq. (]3.9|) taking into account that the coefficient N„ 



(c) 



vanishes. For this particular specialisation Eq. (3.9) contains only terms of the form 

S n (X 1,n ~ 1 \Y 1 ' n ). In this way we use formula (4.12) and set A^ = fii — 7 to obtain the 
relation 



n-l 



n-l 



J V{X^ n - l \Y^ n - 1 ) + J2 K i B)v ( xl ' n ~ 1 \ Y ?' n ~ 1 ) + ^Kf ) V{X° i ' n - l \Y l > n - 1 ) = . 



i=i 



i=l 



(4.15) 



Thus, considering Eq. (4.14) in addition to (4.15), we can see that both functions V and 



W also obeys the Eq. (3.9) under the mappings L \- > L — 1, n ^ n — 1 and /ij /ii+i- 



Step 9. Next we consider Eq. (3.9) with Aq = \i\ — 7 and make use of the relations 



(4.6) and (4.3). By doing so we obtain the expression 



S n {X^ n \Y^ n ) =J2t v ( x l' n \ Y j 1,n ) W(X± n \Y*> n ) , (4.16) 
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where 




Figure 1: Interdependence among the Steps 1-10. 







n 



a(Xf - ^ k )b(n k - Xf) 



b(Xf - fJ,i)b(jii - Xf) a (fjti - fik)a(fi k - Mi) 



x 



n«(^-Mi) 



*;=i 



W-Af)ll 

A a(Af - /x fc )6(/i fc - Af ) 
6(Af - Hi)b(Xf - hi) J- = l a(/ii - fi k )a(fik ~ Mi) 



TT Q(Afc-Ml) TT «(Af-/il) 

6 (Afc - Mi) M^Af-Mi) ' 



J 7 



a(Af - A£) 
Afc) 



(4.17) 



Step 10. For completeness we also set Aq = Mi m Eq. (3.4) and consider the expressions 



(4.10) and (4.12). This procedure yields the same formula (4.16). 
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The interdependence among the Steps 1-10 is schematically depicted in Figure [TJ 
and at this stage we have already gathered the main ingredients required to obtain an 
explicit expression for the function S n . For instance, in the Step 9 we have obtained a 
formula expressing the function S n as a linear combination of auxiliary functions V and 



W defined respectively by (4.2) and (4.5). This process can be thought of as a separation 



of variables induced by the special zeroes described in Lemma ([2]). On the other hand, 
in the Steps 3 and 8 we have shown that the function V satisfy the system of functional 
equations ( 3.4|3.9 ) under the mappings L h \ L — 1, n 4 n- 1 and /ij i-> Mi+i- The 
same holds for the function W as demonstrated in the Steps 4 and 7. Thus, since V 



and W are polynomials of the same order, the linearity of the system of equations (3.4 
3.9) tells us that W = aV where a is a constant. This uniqueness property employed 



in our argumentation has been proven in 27 under very general conditions. Moreover 



since S n is essentially a multivariate polynomial and so is the auxiliary function V due 



to (4.2), the residues of (Oij + a$j t i)/J~ must vanish. This condition tells us that 
a = 1. Thus for n < L the formula ( |4.16[ ) allows us to obtain the function S n up to 
an overall multiplicative constant starting with the solution of (3.4, 3.9) for the case 
n = 1. The solution Si(Af |Af) has been obtained in Appendix [E and in what follows 



we shall demonstrate that the iteration procedure described by (4.16) can be mimicked 
by a multiple contour integral. 



Multiple contour integral. According to the above discussion the expressions (4.16) 



and (4.17) can be rewritten as 



where 



s n (x x > n \Y x > n ) = icJ2 ^.,n.Y,'-" ))- 



Ilk=2 a (/ i i _ Mfc) a (Mfc - Mi) 
1 



a(Af - //i)6(Ap - Mi)o(Af - ^i)6(Af - Mi) 



f] a (\B u)b(u A ^fT «(Af-Ag) A a(Af-Af) 



fc=i 



n«(Af - , k) Hp k - Af ) n n a(Af " A ' ' 



k=i 



ii 6(A? - A?) ^ 6(Af - Af ) 



(4.18) 



Moreover, the formula (4.18|) suggests that the function S n can be expressed as 
5 n (X 1 ' n |y 1 ' n ) = 



-Qdwidwi H(w 1 ,...,w n \w 1 ,...,w n ) 



11 2ivr 2i7r Yll J=1 b(w t - Af )6(^ - Af ) ' 



(4.19) 



(4.20) 
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where the integrals over the set of variables {w^} are performed around the contours 
enclosing solely the poles at to, = \^ . On the other hand, the contours associated with 
the integration over the set {w.j} contain only the poles at w, 



B C 

also assumed to be independent of the variables A • ' in such a way that the only poles 



A^. The function H is 



contributing to the evaluation of (4.20) are due to the zeroes of ni\=i b( w i ~ rf)b(wi ~ 



Now we can integrate the formula (4.20) over the variables w\ and wi, and by 
doing so we generate the summations over the indexes i and j appearing in (4.20) 



This procedure allows us to look for a term by term identification, and the iteration 



procedure described by (4.18) is realised if we are able to exhibit a function H satisfying 
the condition 



n 



H\ 



I wi - 



n b ( x ? - a * ) n w - n ^ - ^ n ^ - x i 

fc=l k=l k=2 k=2 

k+i k^j 

x /C n itj H(w 2 , • • • , w n \w 2 , ...,w n ). (4-21) 



In its turn the function H consists of the function H, up to an overall multiplicative 
constant, under the mappings L \- > L — 1, n \- > n — 1 and /i, (-)■ /ij + i Q Also, it is 
important to remark here that the relation (4.21) needs to be valid for any i,j £ [l,n] 
and that is required to hold only when integrated according to ( 4.20[ ). Thus, under 

these requirements, we simply need to consider the relation (4.21) under the mappings 

c 



Af i->- wi, Xf i->- wi, A^' i— > Wk for k ^ i and Af h- > Wk for k ^ i to obtain the relation 

H(w 2 , . . • ,w n |w 2 , ...,w n 



H(wi, . . .,w n \wi, ...,w n ) 



b{w x - fii)b(wi - /ii) 



] J b(w! - w k ) 2 b(wi - w k ) 2 



k=2 



X 



X 



n 

k=l 



a{w k - /ii, 



Ilfc=2 a (/ i i ~~ / i fc) a (/ i fc - 

L 

n - Mt )6( w - Wl ) n "I" 1 _ 2[ n 



n 



a{w k - /ii, 



b(w k - /ii) b(w k - /ii) 



-i 



a(wi - w k ) A a ( w fe _ w i) 



k=l 



, :2 K^l -W*) ^ &(^fc - Wi) 



]^[ a(w x - Hk)Knk ~ wi) J | 



fe=i 



a(w fc - wi) -A- a(wi - w k ) 
. 2 b(w k - wi) 6(wi - w k ) 



(4.22) 



Now the expression (4.22) can be readily iterated once we know the function H{w\\wi). 
From formula (E.8) we can immediately read that[^] 



HiwAwx) 





l^k)b(w 1 - /i fc ) - 


nt 


1 a(wi - Hk)b{wx - /i fc ) 


b{w\ 


- Hi)b(wx - /ii) 


b(u>i 


-/ii) a(wi-m) 
-Hi) 6(wi— jui) 





(4.23) 



1 Strictly speaking, the relation (4.21) is only required to be valid when integrated as 
■■§[ ] 11?=2 dw i d ™v 

2 Alternative contour integrals expressions are also possible for the case n = 1. Here we have chosen 



(4.23) in order to have a formula compatible with (4.20) 
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and the iteration of (4.22) leave us with the expression, 



H(wi, . . .,w n \wi, ...,w n ) 



[J b(wi - Wj) 2 b(wi - Wj) 2 a(wj - fa)a(wj - fa) 



nr=i K w i - i*i)K w i - & 



i=i 



(4.24) 



where the functions Ri and are given by 

a(w k - fa) j-r a(w k - fa) 



A, 



n 



b(w k - /ij 



n 



- fa 



k=i k—i 

n -p-r a(wi - w k ) a(w k - wi) 
a(u, - ^(w - «*) jl 6N _„, 4) K^.ffl.) 



k=i 



a[w k - Wjj alii;* - WfcJ 



fc=j fc=i+i 



(4.25) 



The formula (4.24) already takes into account the asymptotic behaviour described in 
Lemma HI 



4.2 On-shell formula 



In what follows we shall consider the iteration of the relation (4.18) when the set of 
variables Y 1,n are constrained by the Bethe ansatz equations (2.13). In that case we can 



see from (4.19) that the term will be the most affected and the determination of 
the function H(wi, . . . , w n \wi, . . . , w n ) will require some extra effort. For instance, we 
shall need to promote the function H to H^ s \w S) . . . , w n \w s , . . . , w n ) where the index s 
is introduced in order to keep track of how many iterations we have performed. The case 



s = 1 then yields the function H entering in formula (4.20). 



We then follow the procedure described in Section 4.1 keeping in mind that the 



constraint (2.13) should hold at each level of the iteration process. By doing so we find 



the relation, 



#(«), 



~ l ) L Uk=l b (Ws - fa) nLs+1 a ( W k - ^s)b{w k - W s ) 

b(w s - fi s )b(w s - fa) Y\ k=s+1 a(fa - fa)a(fa - fa) 



uk 



w n \w s , ...,w n ) 



x JJ a(w s - w k )a(w k - fa)b(w k - w s ) 

k=s+l 

L 

a{w, 

k=s k=s+l 
1 L 



n 

,k=s 



a{w k - fa) 
b(w k - fa) 



n 



a{w k - fa) 
b(w k - fa) 



x 



(-1)" 
Ull\ b (w s - fa) 



n 



fa) I [ a{w k 



w. 



s-l 



] J b(w s - fa) ] [ a(w s - w k ) Yl 



k=s 



k=s+l 



k=l 



a{w s - w k ) 
a(w k - w s ) 
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x# (s+1) (> s+1 , . . .,w n \w s+1 , ...,w n ), 



(4.26) 



which needs to be iterated starting from s = l up to s = n — 1. Then at t he la st step 
we shall also need the function H^ n \w n \w n ) which is readily obtained from (E.8) taking 
into account the relation (2.13). Thus we have 

HW(w n \w n ) = {-l) n c 



n 2 YlLlK^n- Hk) 



("If 



0201 



b(w n - ^ n )b(w n - fl r , 
L 

Y\ a(w n - ii k ) 

k=n 
1 



a[w n -fi n ) a[w n -ji n ) 



b(w n - pLn) b(w n -II n )_ 



n-1 



I J b(w n - /i fc ) J ] 



k=n 



a[w n - w k ) 
L a(w k - w n ) 



(4.27) 



and the iteration process above described yields the formula 
H(wi, . . .,w n \w 1 , ...,w n ) = 

nr=i nf=i - n) n^=i - - - w %) 

(_l)(n-l)(i+f) c 2n J _ ./ < 



nr=i 6 (^i - ^ 



]>i 



nr=i - Hi 



OS 



(4.28) 



8=1 



where i?, is given in (4.25) and 



A 



OK 



J^^tUi - /i fc ) JJ a(w fc - 



fc=i+l 



mil a {wi - Vk. 



k=i+l 



i-1 



] J - /ifc) ] [ a(wi - w k ) J J 



L a(u; fc - Wi) 



(4.29) 



The formulas (4.20) and (4.28) constitute an expression for the scalar product S n 



under the constraint (2.13). In that case we have available Slavnov's formula [34] and 



(4.20, 4.28) then corresponds to an alternative representation. Here it is also important 



to remark that it might still be possible to obtain simpler representations for the on-shell 



scalar product from the functional equations (3.4) and (3.9) . For instance, notice that 
the quantity 4>iW + 02 vanishes under the condition (2.13). Thus if we multiply 



Eq. (3.4) by <pi and add it to Eq. (3.9) multiplied by 2 , we are left with the relation 



K S n (X 1 ' n ) + Y,K t S n (X^ n ) = 



(4.30) 



8=1 
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with coefficients Kq = 4>iMq + and Ki 

(4.30) we have omitted the dependence of S T 



+ <f) 2 N^> for % G [l,n\. In Eq. 
with the set of variables Y 1,n since now 



they are completely fixed by the relation (2.13). Moreover, we have checked for small 
values of L and n that Eq. (4.30) solely is able to determine the on-shell scalar product 
up to an overall multiplicative factor. Thus the direct study of (4.30) might still offer 
the possibility of deriving alternative representations. 



5 Concluding remarks 



In this work we have derived functional relations describing the scalar product of Bethe 
vectors for the six-vertex model. The origin of the functional equations is a consistency 
condition between the Yang-Baxter algebra and the highest weight representation theory 
of the 5(2 algebra. More precisely, we have obtained a system formed by two equations 



(3.4, 3.9) whose resolution then yields a multiple integral representation (4.20, 4.24) 



for the aforementioned scalar product. Although a multiple integral formula for this 
same scalar product had been obtained previously in |35| as a specialisation of Baxter's 
solution of the Z-invariant six-vertex model, the structure of the formula obtained here 
is considerably different. 

The resolution of our system of functional equations follows a simple sequence of 
systematic steps which are described in Section 4.1 It is worth remarking that the 



method employed in Section 4.1 consists of an extension of the one developed in 27 ,28 



for the partition function of SOS models with domain wall boundaries. 



Moreover, in Section |4.2| we have also obtained a formula for the scalar product 
(2.12) under the on-shell condition (2.13). In that case it is well known that the six- 
vertex model scalar products are given by Slavnov's determinant formula 34 , and thus 



the formula (4.20, 4.28) constitutes an alternative representation. Furthermore, here 



we have also obtained a single functional equation describing scalar products under the 



on-shell condition (4.30). In this way the direct study of (4.30) might still offer the 



possibility of deriving different representations. 

Recently, there has been a lot of discussion on the possibility of obtaining determinant 
representations for on-shell scalar products in models based on higher rank algebras 
36 39 . In particular, some doubts on that possibility have been put forward in 



Since the method described here is based on the Yang-Baxter algebra, which is the 
common algebraic structure underlying integrable vertex models, it would be interesting 
to investigate the extension of this method for higher rank algebras together with the 



generalisation of the integral formulas (4.20, 4.24 , 4.28). 
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A Polynomial structure 



The structure of the function S n with a given variable Af is completely encoded in the 
operator B(Xf), while its dependence with the variable Xf is described by the operator 
C(Xf). Thus we mainly need to characterise the dependence of the operators -B(A) and 
C(A) with its spectral parameter in order to describe the structure of S n with respect 
to the set of variables {Af } and {Af }. In order to proceed with this analysis it will be 
useful to recast formula ( |2.3 ) as 



r (L) (A) 



A L (X) 
Cl(X) 



B L (X) 



(a.i; 



where the have introduced the index L to emphasise we are considering the ordered 
product of L matrices 7Z a j according to (2.6). In its turn the matrix lZ a j consists of a 



matrix in the space V a whose entries are then matrices acting non-trivially on the j-space 
of the tensor product Vi <8> • • • <8> Yl- Here V,- = C 2 and more precisely we have 



K 



<>j 



7j &j 



where 



a-, 



a 

b 

c 





ft 





c 

b 

a 



(A.2) 



Thus the construction of the operators A, B, C and D described in (2.6) can be imple- 
mented in a recursive manner with the help of the relation 



T {L+l \\) = T (i) (A) TZ aL+1 (\ - fi L+1 ) 



and initial conditions 

A 1 (X) 
Ci(A) 



=«i(A - //i) 

=7i(A - m) 



Bi(X) 
Di(X) 



5i(X- 



Mi) 



(A.3) 



(A.4) 



In terms of its components the recursion relation (A.3) reads 



A L +i(X) = A L (X)a L+1 (X - fi L+1 ) 

B L+1 (X) = A L (X)f3 L+1 (X- fi L+1 ) 

Cl+i(A) = C L (X)a L+1 {X - 

D L+1 (X) = C L (A)/3 L+1 (A-/i L+1 ) 



£l(A)7l+i(A- Hl+i) 
B L (X)5 L+ i(X - hl+i) 
£>l(A)7l+i(A - fiL+t) 
D L (X)5 L+1 (X - /iL+i) 



(A.5) 



which allows us to infer the structure of the operators A(X), B(X), C(A) and -D(A) with 
respect to the spectral parameter A. For that it is important to keep in mind that the 
spectral parameter A only enters those operators through the functions a and b since c 
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depends only on the parameter 7. In terms of the variables x = e and q = e 1 we have 
a(x) = 2~ 1 x~z(xq — q^ 1 ), b(x) = 2~ x x~ 2 {x — 1) and c(x) = 2~ 1 (g — and in this way 
the entries of lZ a j can be written in the form 



Ctj [x ) 



1j[ x ) 



--2- l x-*P a {x) 



/3j(x) 
Sj(x) 



2- l x~^PAx) 



(A.6) 



The terms P a (x) and Ps(x) in (A.6) are polynomials of order 1, while Qp and Q 7 are 
constants. Under these considerations the iteration of (A. 5) tells us that the operators 
P>(\) and C(A) are of the form, 



x 2 Pb{x) 



C(A) 



X 2 P C [X) 



(A.7) 



where Pb(x) and Pc(x) are both polynomials of order L — 1. Thus the definition (2.12) 
combined with the property (D.3) immediately tell us that the scalar product S n is of 
the form 



S n — Y\_( X f X f) 2 Sn( x l \ ■ ■ ■ J X n \ X 1 > ■ ■ ■ > X n ) 



(A.8) 



where re, 



-,2A 



and is a polynomial of order L — 1 in each one of its variables. 



B Special zeroes 



The method considered in Section 4.1 for solving the system of Eqs. (3.4, 3.9) relies on 
the identification of certain zeroes of S n as function of the set of variables {Xf } and {Af }. 
This is due to the fact that we are interested in a polynomial solution, and as such it can 
be characterised by its zeroes. In what follows we shall proceed with the identification of 
those zeroes for particular values of n for illustrative purposes, and subsequently consider 
the general case. 



B.l Case n = 2 

The coefficients N[ and AQ vanish for the specialisation Af = \i\ and Af = fi\ — 7. 
Thus for this particular choice of variables Af 2 > Eq. (3.4) simplifies to 



M Q S*{n x ,m - 7|Af , Af ) + N?S 2 {ji x , ^ - 7 |A , A*) + N^S 2 (^, ^ - 7 |A , Af ) = . 



tB 



tB< 



(B.l) 



Next we set Aq = /ii — 7 noticing that Eq. (B.l) does not simplify significantly. Nev- 



ertheless, the function Mo then consists of a single term. In this way we obtain the 
following relation, 



n[ b) 

-j^S 2 (fM 1 ,fM 1 - 7I//1 - 7, Af ) 
ATf) 



Mn 



(B.2) 
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Eq. (3.9) also simplifies similarly under the specialisation Af = Hi and \ 2 — Hi~ 7. In 



that case it reads, 

M S 2 (hi, m - 7 | Af , Af ) + N b S 2 (hi,Hi - 7IA0, Af ) + N^S^i, Hi - j\Xq, Af ) = 



(B-3) 



and we can substitute (B.2) into ( B.3[ ) to obtain 
QqS 2 (hi, Hi - 7|/ii - 7> A ) 



+<5i-S , 2(a*i, Aii - TlA*i - 7) Ai ) + Q 2 S 2 {hi, Hi - l\Hi - 7> A 2 ) = . (B.4) 

In terms of variables xf = e 2A * s , the functions are rational functions whose explicit 
forms are not enlightening. Nevertheless, they possess very non-trivial zeroes and by 
setting Af = {r^ | Qi{r,j) — 0} for i — 1, 2 we are left with the relation 



Qo\\B =n 

A?=r, 



£2(^1, a*i - j\hi - 7> Ao) = . 



(B.5) 



Thus we can conclude that S 2 (hi, Hi ~ i\hi — 7> Ao) vanishes and by taking this result 
back to (B.2) we obtain S 2 (hi, Hi ~ 7|Af , Af ) = 0. The same procedure can be employed 
considering initially the variables Af instead of Af , and by doing so we end up with the 
property S 2 (\<? , X 2 \hi, Hi~l) = 0. 

B.2 Case n = 3 



We set Af = Hi and Af = /ii - 7 in Eq. (3.4) such that n[ c) = N^ C) = 0. In that case 



(3.4) simplifies to 

M S 3 (hi, Hi ~ 7, Af [Af , Af , Af ) + ivf ^(Aq, hi,Hi~ 7|Af , Af , Af ) 
+n[ b) S 3 (hi, hi ~ 7, Af I Ao, Af , Af ) + N^S^hi, Hi ~ 7, Af |A , Af , Af ) 



+Ni B) S 3 (Hi, Hi ~ 7, Af I Ao, Af , Af ) = 



Cl 



(B.6) 



Next we set Af = Hi ~ 7 an d Af = Hi m (B.6) and obtain the relation 



M S 3 (hi,Hi -7,Af \hi -7,/ii,Af) + Nf i 3s(X , Hi, Hi ~ l\Hi -7,£*i,Af) 
+iV 3 (B) 5' 3 (// 1 , /i! - 7, Af I A , A*i - 7,^i) = . (B.7) 



In addition to that we set Aq = //1 — 7 in (B.7) which then yields the following relation 



M S 3 (hi, Hi ~ 7, Af \hi ~ 7, Hi, Af ) + nPs 3 (hi -^,Hi,Hi~ t!a*i ~ 7, Af ) 
+N^ B) S 3 {fj,i, Hi - 7, Af I /ii - 7, /ii - 7, /ii) = . (B.8) 



Furthermore, the relation (B.8) simplifies to 



S3O1 - 7, Hi, Hi ~ l\Hi ~ 7, Hi, A) = S 3 (hi, Hi ~ 7, Ml 1 ! ~l,Hi- 7, A*i) 



(B.9) 
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for Af 



Af = A since the coefficient M is proportional to 6(Af — Af ). Then we can see 



that Eq. (B.8) combined with the property (B.9) yields the relation 



Sa{m, Mi - 7, Af|/ii - 7, Mi, A3 ) 



A'. 



M 



-S^Mi - 7, Mi, Mi - 7lMi - 7, Mi, Af ) 



AT, 



(B) 



TT- S 3 (m -7,Mi,Mi -7lMi -7,Mi,Af) 

iW 



(B.fO) 



which can be substituted back into (B.7) considering the symmetry property (|3j). We 
are then left with an equation of the form 



QoS^Mi - 7, Mi, Mi - 7lMi - 7, Mi, Aq) + Q 3 ^(Mi ~ 7, Mi, Mi - 7lMi - 7, Mi, A 



+Q3*5 , 3(mi - 7, Mi, Mi 



7I/Z1 — 7, Mi, A^) = 



(b.ii; 



where Qq and Q 3 ,Ll are rati onal f unctions in the variables xf ,c . Then we can set Af 



c 



B,C fr B,C 



0} in ( |B.ll ) and conclude that Ss(fii — 7, [ix, fix — 7|mi — 7, Mi, Ac 



0. Considering (B.10), the latter also implies in 

^(Mi, Mi - 7, A^Imi - 7, Mi, A 







(B.12) 



Next we set Ao = [ii and Af = ji\ — 7 in Eq. (3.4) considering the property (B.12). By 
doing so we obtain the following relation, 



M S 3 {vi, /ii - 7, Af|/ii - 7, Af , Af ) + A^ 3 l ° J 5 3 (/ii, /ii, /ii - 7I//1 - 7, Af , Af ) = . 



(CO, 



(B.13) 



The coefficient M in (B.13) vanishes for A3 7 



/ii — 7 while is finite. Thus for 

this particular specialisation of Af we obtain S 3 ([ix, fix, Mi — 7lMi — 7, Af , Af ) = 0, and 
consequently we can conclude that 

^(Mi, Mi - 7, Af |mi - 7, Af , A 



. 



(B.14) 



So far we have only considered Eq. (3.4) and the next step is to set Af = fix, 
\% = fi\ — 7 and Ao = /ii — 7 in Eq. (3.9). We also take into account the property (B.14) 
and by doing so we obtain the relation 



M £ 3 (Mi, fix - 7, Af I Af , Af , Af ) + Nt'Ssitix - 7, Mi, Mi - 7|Af , Af , Af ) = . (B.15) 



HQ 



B \B \B\ 



The function M possesses non-trivial zeroes and we can set Af in (B.15) such that M 
vanishes in order to conclude that S 3 (/Jx — 7, Mi, Mi — 7|Af , Af , Af ) = 0. Then from 
(B.15) the latter implies in 



> S 3 (Mi,Mi-7,A 3 7 |Af,Af,Af) = 0. 



(B.16) 



The procedure above described can also be performed with the specialisations of the 
variables Af and Af exchanged allowing us to conclude that 



^(Af,Af,Af| yUl ,/i 1 - 7 ,Af)=0 



(B.17) 
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B.3 General case 



B,C 



\f_^ + 7 for j G [3, n — 1] and collect the results at each step. At the final step 

fix — j. This 



We consider Eq. (3.4) under the specialisations X\ 



c 



Af 



A*l> Af 



Af = Hi — 7 and 



A 

the only non- vanishing coefficients are M and Nn ' , and we then set Aq 
procedure then yields the formula, 



S n (X\Y) 



N., 



(C) 



N, 



(B) 



, r -5 n (x*|F) - -^s n (x|y*) 



(B.18) 



where 



y 

Y* 



{/ii, /ii - 7, Af + (n - 3)7, . . . , Af + (n - jh, • • • , Af} 
{//i - 7, /xi, Af + (n - 3)7, . . . , Af + (n - Af } 

{/ii - 7, /ii, /ii - 7, Af + (n - 3)7, . . . , Af + (n - j)7, . 
{/ii - 7, n x - 7, /ii, Af + (n - 3)7, . . . , Af + (n - j)7, • 



An +7} 
Af + 7} 



(B.19) 



The relation (B.18) can now be substituted back into the previous steps leading to it, 



and an analysis similar to the one employed for the cases n = 2, 3 allows us to conclude 
that 



s n (m,ni -7, a 



c 

3 ) 



Af|/i 1 - 7 ,A^,...,Af) = 0. 



^B.20) 



Next we set Af = fii, Af = /ii — 7 and Aq = /xi — 7 in Eq. (|3.9[) taking into account the 



property (B.20). By doing so we obtain the more general condition 
S n (ni,Hi - 7, Af , . . . , Af |Af , . . . , Af ) = . 



(B.21) 



The procedure above described can also be performed with the specialisations of the 
variables Af and Af exchanged. In that case we then obtain the vanishing condition 



5 n (Af , . . . , Af jj,x - 7, Af , . . . , Af ) = 



[B.22) 



C S n as a doubly symmetric function 



The function S^Af , . . . , Af |Af , . . . , Af ) defined by (2.12) is expected to be invariant 



under the exchange of variables Af <H- Af and Af <H- Af separately. This property is due 



to the definition (2.12) and the commutation relations 



[5(A 1 ),B(A 2 )] = [C(Ai),C(A 2 )]=0 



(CI) 



described in (3.2) and (3.3). Nevertheless, once we assume S n is determined by the Eqs. 



(3.4) and (3.9), it would be desirable this symmetry to be an inherent property of their 



solutions. This is indeed the case and it can be demonstrated as follows. 
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We firstly notice that the coefficients Mq and N- defined in (3.5) are the only 
coefficients exhibiting poles when Ao — > Af . Moreover, those coefficients also satisfy the 
property 



lim M b(X - A 

Ao^Af 



lim Nf b(X - Af ) 

Ao^A? 



cf[a(Xf - H )f[ aiXf Af) 



(C.2) 



Next we integrate Eq. (3.4) over the variable Ao and around the contour Cf enclosing 
solely the variable Af . Thus considering the property (C.2) we obtain the identity 

L n a(X? — A?) 

c J|a(Af - fij) Y\ ^ B _ ^ [^n( A i, • • • > Af |Af , . . . , Af ) 



q (\C \^\\ B \ B \ B \ B \ B 

"JnV A l J • • • i A n \ A i i A l > • • • > A i-1) ' ' ' ' A n 



. 

(C.3) 



The formula (C.3) allows us to conclude that 

S n ( A i , ■ ■ • , A n |Af , Af , . . . , Af_ 1; A^_ 1; . . . , Af ) = 5 n (Af , . . . , Af |Af , . . . , Af ) 



for i G [1, L], which implies the symmetry relation 



(C.4) 



S n (X 1 , . . . , X n I . . . , A i , . . . , Xj , . . . ) — S n (X 1 , . . . , X n | . . . , Xj , . . . , Aj , . . . ) . (C.5) 

Now in order to demonstrate an equivalent symmetry relation with respect to the set 
of variables {Af }, we integrate Eq. (3.4) over the variable Ao and around the contour 
Cf containing solely the variable Af . We also take into account that the coefficients Mq 
and are the only ones possessing poles when Ao — > Af . In this way we obtain the 
identity 



c Y]_ a ( A ? - N) \\ h f\c _ \c\ [^n(Af) ■■■■> Af |Af , . . . , Af ) 

5=1 



;iKAf-Af) 

q (\C \C \C \C \C\\B 

°n\ A i ) A l j • • • j A i-li A i+li ■ ■ ■ j A n \ A \ i 



Af )] = 



(C.6) 



with the help of the property 

lim M 6(A - A 

Ao^Af 



- lim NP b(X - Af) 

Ao-^Af 

-cf[a(Xf - fi^fl f^ 

j=l 3=1 ^ J * ' 



(C.7) 
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The property S n (Xf , Af , . . . , Af_ l5 Af +1 , . . . , A?|Af , . . . , A* ) = S n {\? , 
is then valid for i £ [1 , L] which allows us to conclude that 



A n |Af 



A*) 



A 



c 



Aj , . . . \A X , . . . , A n ) 



S a 



Aj , . . . , Aj , 



• •• K a;/), (c.8) 



It is important to remark here that the symmetry relations (C.5) and (C.8) have been 
obtained solely from the examination of the Eq. (3.4). Alternatively, we could also have 
derived the same symmetry relations from a similar analysis of the Eq. (3.9). 



D Asymptotic behaviour 



Considering the variables x = e , q = e 7 and the conventions described in (A. 2) and 



(|A.2|), in the limit x — > oo we find 

K 



a 



i i 
X2q2 



q-q v h 

7 ~ X 



2 

i i 



K- 1 



where the operators K and X ± are explicitly given by 



K 



q? 
q~ 



X z 



1 ±1 
2V1T1 



In this particular limit the relations (A. 5) can be easily iterated and we obtain 

L 

k 



B(X) 
C(A) 



(q-q x ) h=L 

- -q 2 x 2 e — 



— — -q 2 x 2 e ^=i^ 

2 y 



fc=i 



fc=i 



with operators Pp being defined as 



P ± 



3-1 L 

,Jc - ® X ± ® (g) i^ 1 . 

k=l k=j+l 

In their turn the operators K ±x and X 1 * 1 satisfy the g-deformed su(2) algebra, 

KX ± K~ 1 = q ±1 X ± 

[x\x-] = 



K 2 -K~ 



q-q 1 



which allows us to demonstrate the relations 



pa pb _ (a,b) pb pa 
i 3 q 3 i 

(p-) 2 = o 



I* < 3 j 



(d.i; 



(D.2) 



(D.3) 



(D.4) 



(D.5) 



(D.6) 
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with symbols (a, b) being denned as 

(a, b) = 



2 a = ±,b = ± 
-2 a = ±,b = T 



(D.7) 



Now with the help of the relations (D.6) and considering Xj = e 2A % in the limit 
Xi — )■ oo we obtain the expression 



ri B ' 



and similarly 



-l\n 



2 nL 



l( ^T 11 - n ( n - 1 )fr).!] 9 2'=~ n ^'- i /,f 



-if ~ ■",//!),,■( ] J., 



2 



x E 

l<ai<---<a n <L 



En 



n 



i<i<™ 



(D.8) 



i=i 



2 nL 



n 

t=l 



.1' ■ 



E 



En 
3=1 Ma i 



l<ai<---<a n <L 



n p i- 

l<?<n 



(D.9) 



The term [n\] q 2 appearing in (D.8) and (D.9) corresponds to the g-factorial function and 
it is defined as [n!] ?2 = 1(1 + qrfl + q 2 + 7f) . . . (1 + q 2 + ■ ■ ■ + g 2 ^" 1 )). 

The relations ( D.8[ ) and (D.9) can now be combined according to (2.12) to obtain the 
following formula, 



/ — l\2n n 

S n (x°, ■ ■ • ,XZ\X° , • • .,0$) ~ ii^^l)-^!)^^^ 11^?)^ 



i=l 



e e^^+^( | n p a + n p-| ) 



l<ai <■ ■ <a n <L 

l<6i<-<6 n <i 



(D.10) 



in the limit sf'* 7 oo. The next step is to compute the quantity (0 1 P+ ff P 6 . |0) 



appearing in (D.10), and for that we notice that under the constraint 1 < a\ < ■ ■ ■ < 
a n < L and 1 < b\ < ■ ■ ■ < b n < L the aforementioned quantity will contribute only 
when bj = aj. In that case we have (0| P+ P~ |0) = g n ( n_1 ) and the relation 



l<J<n 3 l<j<n 



(D.10) then simplifies to 



S n ~ 



(Q Q »(£-») r_ |12 _-2n£tiW 



22nL 



-5 ' ' ™La e 



E 



l<ai<---<a n <L 



(d.ii; 
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E Off-shell solution for the case n = 1 



The relation (4.18) offers a way of building the solution of (3.4, 3.9) recursively. More 



precisely, the formula (4.18) establish a relation between the scalar product S n and an 
auxiliary function V consisting essentially of S* n _i for a lattice of length L — 1. Thus, 
since we are considering n < L, at the final step of this iteration procedure we shall 
need the solution of Eqs. (3.4) and (3.9) for the case n — 1. In that case our system of 



functional equations explicitly read 

M Si(A? |Af ) + Ni B) S 1 (X^\X ) + TVf^rOolAf ) 
M £i(A? |Af ) + N^StiX? |A ) + NPs 1 (X \\f) 

with coefficients 

c(A - Af ) 6(Af - A?) A _^ at(b) 








(E.l) 



M 



M 



Af) 

'fe(Ao-Af) 6(A -Af) 

c(A - Af ) 6(Af - Af ) 
'fe(Ao-Af) 6(Ao-A^ 

Af x L 



c(A c 



6(A - Af ) 



Y[a(X - H ) Nl 

3=1 
L 

H 6(a -m) -v; 

3=1 



c(Af 



OB) 



(C) 



6(Af - A ) 

c(A - Af ) 
fc(A - Af ) 

c(Af - Ap) 
6(Af - A ) 



Ao) n^ 

3=1 



B 



Pi, 



flj) 



3=1 

L 



HKxf-^) 



3=1 



(E.2) 



and the solution of (E.l) can be obtained as follows. By eliminating the term Si(Ao|Af ) 
from (E.l), we obtain an equation involving only the terms Si(Af|Af) and Si(Af|A ). 
More precisely, we are left with the identity 



(E.3) 



6(Af - 


-Af^Af 


Af) 






6(Ao- 


-Af)5i(Af 


Ao) 




rr L a( ~ x ?- 


-Hj) y\L 
-N) lL 3= 1 


6(Af 
6(A? 


-Mi) 
~Mi). 




■p-j-L a(A - 
|llj=l a (Af- 


~Mi) TT L 
-Mi) Ali=l 


fe(A - 
6(Af- 


-Mi) 
-Mi). 



where the variable Af now assumes the role of a parameter. The relation (E.3) is then 
readily solved by 

Si(Af|Af) 



F(Af) 



n 



n 



6(Af - /*) 



L a(Af - fij) Al 6 ( A f _ 



(E.4) 



where F is an arbitrary function. Next we substitute the expression (E.4) into the first 



Eq. of (|E.l|), and after eliminating an overall factor we are left with the relation 

Hj)b(X - Hj) A a(A - fij)b(\f - fij) 



F(Af; 

-F(Ao) 



n 

.3=1 



o(Af 



B 



6(Af - 



n 



6(Af - /ij) 



n 



a(Af — fij)a(Xf — fij) T-r a(Af — fij)b(Xf — /ij) 



a(A - fij) 



n 



6(A - /ij) 



0. (E.5) 
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As far as the Eq. (E.5) is concerned, the variable Af is merely a parameter which can 
be set at our convenience. In particular, for the specialisation Af = — 7 the Eq. (E.5) 
simplifies to 



^(A c 



F{Xi 



n,=i a (^o - ^)&(A - fij) n,=i a ( x ? - - Hj) 



(E.6) 



Eq. (E.6) can then be easily solved and we find 

L 



F(X) = Q Y[ a(X - Hj)b(X - ^) 
3=1 



(E.7) 



where Q is a A independent parameter. Then the combination of (E.4) and (E.7) yields 
the expression 



^(AflAf 



c(Af - Af ) 
6(Af - Af ) 



' L L 

H a(Af - nj)b(\? - m) - J] a(Af - H )b(Xf - N ) 

.3=1 3=1 



(E.8) 



where the constant i? has been fixed by the asymptotic behaviour (D.ll). 
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